Yukawa sector in non-supersymmetric renormalizable SO (10) 
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We discuss the ordinary, non-supersymmetric SO(10) as a theory of fermion masses and mixings. 
We construct two minimal versions of the Yukawa sector based on 126h and either 10h or 120h. 
The latter case is of particular interest since it connects the absolute neutrino mass scale with the 
size of the atmospheric mixing angle 9a ■ It also relates the smallness of V c b with the largeness of 6a ■ 
These results are based on the analytic study of the second and third generations. Furthermore, we 
discuss the structure of the light Higgs and the role of the Peccei-Quinn symmetry for dark matter 
and the predictivity of the theory. 



I. INTRODUCTION 

SO(10) grand unified theory p| is probably the best 
motivated candidate for the unification of the strong 
and electro-weak interactions. It unifies the family of 
fermions; it includes the SU(4)c quark-lepton symme- 
try and the left-right (LR) symmetry Q in the form 
of charge conjugation as a finite gauge symmetry; it pre- 
dicts the existence of right-handed neutrinos and through 
the see-saw mechanism Q| offers an appealing explana- 
tion for the smallness of neutrino masses. Due to the 
success of supersymmetric unification, and the use of su- 
persymmetry in controlling the gauge hierarchy, most of 
the attention in recent years has focused on the super- 
symmetric version of SO (10). However, supersymmetry 
may not be there. After all, it controls the Higgs gauge 
hierarchy, but not the cosmological constant. The long 
standing failure of understanding the smallness of the 
cosmological constant suggests that the unwelcome fine- 
tuning may be necessary. Our fine-tuned world can be 
viewed, in the landscape picture simply as a selection 
criterion among the large number of degenerate string 
vacua. Or it could be that the cosmological evolution of 
the universe selects a light Higgs doublet [jfl. If so, the 
main motivation behind the low-energy supersymmetry 
would be gone. However one could still hope that the 
issue might be settled by the internal consistency and 
the predictions of a well defined grand unified theory. In 
other words, it is possible that simplicity and minimality 
at the GUT scale requires a specific low-energy theory. 

What about grand unification without supersymme- 
try? At first glance, one may worry about the unification 
of gauge couplings in this case. Certainly, in the minimal 
SU(5) theory, the gauge couplings do not unify without 
low-energy supersymmetry. What happens is the follow- 
ing: the colour and weak gauge couplings meet at around 
10 16 GeV, an ideal scale from the point of view of the pro- 
ton stability and perturbativity (i.e., sufficiently below 
M Planck)- The problem is the U(l) coupling. Without 
supersymmetry it meets the SU(2) r, coupling at around 
10 13 GeV Jj|; with low-energy supersymmetry the one- 
step unification works as is well known . 



On the other hand, the fact that neutrinos are massive 
indicates strongly that SU(5) is not the right grand uni- 
fied theory: it simply requires too many disconnected pa- 
rameters in the Yukawa sector 

HIT The SO(10) theory 
is favored by the neutrino oscillation data. Most inter- 
estingly, SO(10) needs no supersymmetry for a successful 
unification of gauge couplings. On the contrary, the fail- 
ure of ordinary SU(5) tells us that in the absence of super- 
symmetry there is necessarily an intermediate scale such 
as the left- right symmetry breaking scale Mr. Namely, 
in this case the SU(2)l and SU(3)c couplings run as 
in the standard model or with a tiny change depending 
whether or not there are additional Higgs multiplets at 
Mr (recall that the Higgs contribution to the running is 
small). However, the U(l) coupling is strongly affected 
by the embedding in SU(2)^ above Mr. The large con- 
tributions of the right-handed gauge bosons makes the 
U(l) coupling increase much slower and helps it meet 
the other two couplings at the same point. The scale 
Mr typically lies between 10 10 GeV and 10 14 GeV (see 
for example 0, [n| and references therein), which fits 
very nicely with the see-saw mechanism. Now, having 
no supersymmetry implies the loss of a dark matter can- 
didate. One may be even willing to introduce an addi- 
tional symmetry in order to achieve this. In this case it 
should be stressed that SO (10) provides a framework for 
the axionic dark matter: all one needs is a Peccei-Quinn 
fl2| symmetry U(l)pg which simultaneously solves the 
strong CP problem. 

This seems to us more than sufficient motivation to 
carefully study ordinary non-supersymmetric SO(10). 
What is missing in this program is the construction of a 
well defined predictive theory with the realistic fermionic 
spectrum. This paper is devoted precisely to this task. 

In particular, the search of the minimal realistic 
Yukawa sector is a burning question. In the absence 
of higher dimensional operators at least two Higgs mul- 
tiplets with the corresponding Yukawa matrices are 
needed, otherwise there would be no mixings. The 
Yukawa Higgs sector can contain 10h, 120h and 126h 
representations, since 
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16x16 = 10 + 120 + 126. (1) 

One version of the theory with only 10h and 126h was 
studied in great detail in the case of low-energy super- 
symmetry [13, Il4l Il5j . In spite of having a small number 
of parameters it seems to be consistent with all the data 
111 Elll El E3 For the type II seesaw it predicts 
furthermore the 1 — 3 leptonic mixing angle not far from 
its experimental limit: \U e 3\ > 0.1 [lSl l2fll| and it offers 
an interesting connection between b — r unification and 
the large atmospheric mixing angle 0, l^]. Last but 
not least, it also predicts exact R- parity at low energies 
[23l |24| leading to the LSP as a candidate for the dark 
matter. 

Thus, a first obvious possibility in ordinary SO (10) 
is to address the model with 10h + 126h, and to see 
whether or not it can continue to be realistic. 

There is another interesting alternative: 120h and 
126h- We limit ourselves to the analytic study, which 
requires ignoring the effects of the first generation. The 
main result of this model is the correlation of neutrino 
masses with the value of the atmospheric mixing angle, 
true both in Type I and Type II seesaw mechanisms, or 
in a general case when both contribute to the neutrino 
masses. Furthermore the large atmospheric mixing angle 
fits naturally with the small V c b mixing. 

Both cases require complexifying the Higgs fields 10h 
or 120 h- This in turn calls for a PQ symmetry, in any 
case useful for a dark matter. 

At first glance it seems that the Yukawa sector in a 
non-supersymmetric theory does not differ from the su- 
persymmetric version and thus there is nothing new to 
say. There are however several subtle differences: 1) the 
running of the couplings (gauge and Yukawa) is changed, 
so are the inputs for a numerical evaluation at Mqut] 
2) there are necessarily intermediate scales; 3) if no new 
symmetries are invoked, all the SO(10) representations 
that are not complex like 16 or 126 are real; 4) there are 
radiative corrections to the Yukawa sector, that should 
in principle be taken into account. All these points will 
be discussed below. 



II. THE MINIMAL YUKAWA SECTOR 

In this work, we stick to the renormalizable ver- 
sion of the see-saw mechanism (for alternatives using a 
radiatively- induced see-saw, seejl^j), which makes the 
representation 126h indispensable, since it breaks the 
SU(2)# group and gives a see-saw neutrino mass. By it- 
self it gives no fermionic mixing, so it does not suffice. 
The realistic fermionic spectrum requires adding cither 
10h or 120h- As promised in the Introduction, we will 
go carefully through both possibilities. 

Before starting out, it is convenient to decompose the 
Higgs fields under the SU(2) L x SV(2) R x SU(4) C Pati- 
Salam (PS) group: 



10 = (2, 2,1) + (1,1, 6) 
126 = (1,3, 10) + (3, 1,10) + (2, 2, 15) + (1,1, 6) (2) 
120 = (1,3, 6) + (3, 1,6) + (2, 2, 15) + (2, 2,1) + (1,1, 20) 

As is well known, the 126h provides mass terms for right- 
handed and left-handed neutrinos: 



M VR = (1, 3, 10) Km, M„ L = (3, 1, 10) Yi 26 (3) 
which means that one has both type I and type II seesaw: 

M N = -M VD M~^M vn + M VL (4) 

In the type I case it is the large vev of (1,3,10) that 
provides the masses of right-handed neutrino whereas in 
the type II case, the left-handed triplet provides directly 
light neutrino masses through a small vev j2|| [27| ■ The 
disentangling of the two contributions is in general hard. 



A. Model I: 126 H + 10 H 

In this case the most general Yukawa interaction is 
(schematically) 



C Y = 16 F (10 H Y 10 + 126 H r 126 ) 16 F + h.c. . (5) 

where Y±o and Y\2q are symmetric matrices in the gener- 
ation space. With this one obtains relations for the Dirac 
fermion masses 

M D = M 1 + M , M u = c 1 M 1 + c M , 
M E = -2,Mi + M , M VD = Sc^h + c M , (6) 

where we have defined 



M l = (2,2,15)f 26 F 126 , 

M = <2,2,l)? o r 10 , (7) 

and 



_ {2,2,l)l- _ (2,2,15)? 26 
(2,2,l)f ' Cl (2,2,15)f 26 - W 

These equations, together with © and are the start- 
ing point for the analysis [28} of the fermion spectrum. 
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1. Who is the light Higgs? 

With the minimal fine-tuning the light Higgs is in gen- 
eral a mixture of, among others, (2,2,1) of 10h and 
(2,2,15) of 126h- This happens at least due to the large 
(1,3,10) vev in the term (T26 H ) 2 T26 H f 10 H 

In any case, their mixings require the breaking of 
SU(4)c symmetry at a scale Mps, and it is thus con- 
trolled by the ratio Mps/M, where M corresponds to 
the mass of the heavy doublets. Thus, if M ~ Mqut, 
and Mps <C Mgut> this would not work; we come to 
the conclusion that one needs to tune-down somewhat 
M. Tuning down of the (2, 2, 1) mass cannot have much 
impact on the unification constraints, but (2, 2, 15) is a 
large field and could in principle cause trouble. However, 
its contribution is actually tiny, since the differences in 
the corresponding /3-function coefficients (62 — &3)(2,2,i5) 
and (61 — &2)(2,2,i5) ( m the usual notation) are very small. 
It should nevertheless be taken into account when study- 
ing unification constraints. 

2. The simplest version: real IOh 

If 10h is real, then there is just one SU(2)^ doublet 
in (2,2,1) and thus |(2,2,1)« | = |(2,2,1)? |, namely 
I cq I = 1. The parameter space is thus smaller. Here 
we show that, in the two generation (second and third) 
case with real parameters, such a constraint leads to a 
contradiction with the data. In the physically sensible 
approximation 9 q = V cb = we find 

m c (m T - m b ) - m t (m^ - m s ) m t 

c — « , 

m s m T — m^mb vrib 

clearly much bigger than 1. 

This conclusion is subject to the uncertainties of the 
full three-generation case. Although strictly speaking 
this simple model cannot be ruled out yet, there is an 
indication that a more complicated scenario should be 
considered. 



3. The next step: complex IOh 

If the model with real IOh does fail eventually, one 
could simply complexify it. This of course introduces 
new Yukawa couplings which makes the theory less pre- 
dictive. Certain predictions may remain, though, such 
as the automatic connection between b — r unification 
and large atmospheric mixing angle in the type II see- 
saw. This is true independently of the number of 10 
dimensional Higgs representations, since IOh cannot dis- 
tinguish down quarks from charged leptons. From M VL oc 
^126 , one has M VL oc Mp, — Me- 

It is a simple exercise to establish the above men- 
tioned connection between \rrib\ ~ \m T \ and large 9 a tm] 



for details see 0, . This fact has inspired the care- 
ful study of the analogous supersymmetric version where 
m T ~ mb at the GUT scale works rather well. In the 
non-supersymmetric theory, b — r unification fails badly, 
m T ~ Irrib 29] . The realistic theory will require a Type 
I seesaw, or an admixture of both possibilities. 

4- Axions and the dark matter of the universe 

A complex IOh means, as we said, an extra set 
of Yukawa couplings. At the same time this non- 
supersymmetric theory cannot account for the dark mat- 
ter of the universe, since there are no cosmologically sta- 
ble neutral particles and, as is well known, light neutrinos 
cannot too. It is then rather suggestive to profit from the 
complex IOh and impose the U(1)pq Peccei-Quinn sym- 
metry: 

16 -> e M 16 , 10 -> e- 2la 10 , 126 -» e" 2M 126 , (10) 

with all other fields neutral. The Yukawa structure has 
the form JSJ) with IOh now complex. This resolves the 
inconsistency in fermion masses and mixings discussed 
above, a nd g ives the axion as a dark matter candidate as 
a bonus • 

The neutrality of the other Higgs fields under U(1)pq 
emerges from the requirement of minimality of the Higgs 
sector that we wish to stick to. Namely, 126h is a com- 
plex representation and IOh had to be complexified in 
order to achieve realistic fermion mass matrices and to 
have U(1)pq. It is desirable that the U(1)pq be broken 
by a nonzero (126h), i.e. the scale of SU(2)p breaking 
and right-handed neutrino masses [3l]], otherwise IOh 
would do it an give Mpq w Myy, which is ruled out 
by experiment. Actually, astrophysical and cosmological 
limits prefer Mpq in the window 10 10 — 10 13 GeV |32j . 

Now, a single 126h just trades the original Peccei- 
Quinn charge for a linear combination of U(1)pq, T^p 
and B — L 0, ■ Thus in order to break this combina- 
tion and provide the Goldstone boson an additional Higgs 
multiplet is needed. One choice is to add another 126h 
and decouple it from fermions, since it must necessarily 
have a different PQ charge jSj. An alternative is to use 
a (complex) GUT scale Higgs as considered for SU(5) 
by 34], with Mpq ~ Mqut- which however implies too 
much dark matter or some amount of fine-tuning. 

Of course, the Peccei-Quinn symmetry does more than 
just providing the dark matter candidate: it solves the 
strong CP problem and predicts the vanishing 6. The 
reader may object to worrying about the strong CP and 
not the Higgs mass hierarchy problem; after all, they are 
both problems of fine-tuning. Actually, the strong CP 
problem is not even a problem in the standard model, at 
least not in the technical sense j^]. Namely, although 
divergent, in the standard model 9 is much smaller than 
the experimental limit: 9 <C 10~ 10 for any reasonable 
value of the cutoff A, e.g. 9 « 10~ 19 for A = Mpi anc k- 
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The physical question is really the value of 6. Peccei- 
Quinn symmetry fixes this arbitrary parameter of the 
SM. The situation with supersymmetry and the Higgs 
mass is opposite. Low energy SUSY helps keep Higgs 
mass small in perturbation theory, but fails completely 
in predicting it. If we do not worry about the natural- 
ness we can do without supersymmetry. On the other 
hand, if we wish to predict the electron dipole moment 
of the neutron, U(1)pq is a must, unless we employ the 
spontaneous breaking of P or CP in order to control 6 

IMEllaa. 



introducing a U(1)pq global symmetry, which provides 
as a byproduct a dark matter candidate. 

The type I seesaw contribution due to right-handed 
neutrinos gives the light neutrino mass matrix 



Mpj = -Ml D M-^M VD oc 9c?Mi - c\M 2 M^M 2 (13) 
whereas the type II contribution reads 



Mtf oc Mi. 



(14) 



B. Model II: 126 H + 120 H 

Instead of 10h one could use 120h- Since Y\2Q is 
antisymmetric, this means only 3 new complex couplings 
on top of Yi26- 

120h as an addition to 10h was partially studied some 
time ago, for charged fermions only [l^. Recently it 
was readdressed J25J in the context of the radiative see- 
saw mechanism |39(. The analytic study for the second 
and third generation gives b — r unification, and small 
quark and large leptonic mixing angle. We recall that 
this model is tailor-fit for the strongly split supersymme- 
try |4£j with light gauginos and higgsinos and very heavy 
scalars [4l| . 

The analysis here is quite similar to the case of 10h 
and 120h due to the fact that 10h and 126h have 
symmetric Yukawas. There is an important difference 
though, since (2,2, l)i is traded for (2,2, 15)i26- 

The Dirac mass matrices at the grand unification scale 
take the form 

The Dirac mass matrices at the grand unification scale 
take the form 



Md = Mi +M 2 
Me = -3Mi + c 3 M 2 



M v = c 1 M 1 +c 2 M 2 , (11) 
M VD = -3ciMi +C4M2 



where Mi and c\ are defined in J7J|,(|SJ), and: 

M 2 = yi2o((2,2,l)f ao +(2,2,15)f 20 ) , 
<2,2,1)? 20 + <2,2,15)? 20 



('2 



C3 



c 4 = 



(2,2,l)f 20 + (2,2,15)f 20 ' 
(2,2,l)f 20 -3(2,2,15)f 20 
(2,2,l)f 20 + (2,2,15)f 20 
(2, 2, 1)^0-3(2,2,15)^ 
(2,2,l)f 20 + (2,2,15)f 20 



(12) 



The case of real 120# reduces to similar constraints al- 
ready encountered in the real case of model I (see also 
|42j for an analysis using SU(5) decomposition). For real 
bidoublets the definitions (|12fl constrain all three c,; to 
a same order of magnitude. But this contradicts the re- 
quirements for small second generation masses of charged 
leptons (03 ~ 3) and of up quarks (c 2 « m t /rat>). In other 
words, similarly as in model I, there is a need to com- 
plexify the Higgs fields. This is again best achieved by 



1. Two generations case: analysis and predictions 

In spite of only two Yukawa matrices, the above system 
of equations is rather complicated and requires painstak- 
ing numerical studies. Before plunging in this computa- 
tional tedium, it is certainly useful if not indispensable 
to get a physical insight through analytical arguments. 
Following the successful approach that was adopted by 
us before (for the case of a supersymmetric SO(10) the- 
ory with 10h and 126h 0, E2) we focus here on the 
study of the second and third generations, with the nat- 
ural expectation that the effects of the first generation 
can be treated as a perturbation. 

In the basis where Mi is diagonal, real and non- 
negative: 



, , , sin 2 6 

Ml oc _ 2 
V cos^i 



(15) 



the most general charged fcrmion matrix can be written 
as: 



M f = [i f 



sin i (sin 6* cos # + e/) 

-i(sin0cos0 + e/) cos 2 



(16) 

where / = D,U,E stands for charged fermions and e/ 
vanishes for negligible second generation masses. In other 
words |e/| oc m 2 /m| (see below). Furthermore the real 
parameter [i f sets the third generation mass scale, made 
explicit below. 

We determine next the unitary matrices Lf and i?/, 
which diagonalize Mf in the physically relevant approx- 
imation of small |e/|. We obtain 

M = fl}-diag{-/i/e/ sia20,/i/(l+e/sm20)}-£/+0(|e}|) 

(17) 

where 



L f = 



—i cos 29 €f 



-i cos 28 ef 
1 



cosf 
-i sin( 



—1 smt 
cos 



(18) 
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and Rf is given by the same expression with i — * —i. 
From eq. (|17f) . at the leading order in |e| we get: 



Iwl 

sin20|e/| 



m 



3 ! 



(19) 
(20) 



The phases of the three e/ parameters are not deter- 
mined, while the meaning of the angle 9 will be clear in 
a moment. Now we are in the position to state the three 
predictions of this theory regarding 1) neutrino masses, 
2) the relation between bottom and tau masses, and 3) 
the quark mixing angle V c b- 

1) We begin with the predictions regarding neutrino 
masses. By using eqs. (|13|) and 11141 and an explicit form 
of the 2x2 matrices one concludes that type I and type 
II seesaw lead to the same structure 



M* N 



MU cx Mi 



(21) 



In the selected basis the neutrino mass matrix is diagonal. 
We see that the angle 9 has to be identified with the 
lcptonic (atmospheric) mixing angle 9a up to terms of 
the order of |ejs| ~ m^/rrir. For the neutrino masses we 
obtain from i|15|) 



cos 29 a 



1-sva 2 29 A /2 



0(\e\ 



(22) 



This equation points to an intriguing correlation: the 
degeneracy of neutrino masses is measured by the max- 
imality of the atmospheric mixing angle. What about 
numerical predictions? Clearly, without including the ef- 
fects of the first generation and the impact of the running 
from the GUT to the weak scale, no precise determina- 
tion can be made. It may be illustrative though to give 
an estimate in case this formula were to remain approx- 
imately valid. The value of 777-2 could not then be too 
small: e.g., with the value Am A — |m§ — m|| « 2.5- 10 -3 
eV 2 and the 99 % CL limit 9 A = 45° ± 9° from Q one 
would get m 2 > 30 meV. On the other hand, there is an 
upper limit from cosmology and neutrinoless double beta 
decay, which (depending on the selected data and analy- 
sis) varies from 0.14 eV to 0.5 eV, see again ^3j. Clearly, 
the larger the limit, the closer one can be to 9a = 45°. 
The analysis in 3) below suggests though that 9a should 
be as far as possible from the maximal value, i.e. that 
neutrinos should be as hierarchical as allowed by l|22|l . 

2) The second prediction regards the ratio of tau and 
bottom mass at the GUT scale: 



777 T 



3 + 3 sin 20a Re[e E - e D ] + 0(\e 2 \) 



(23) 



At first glance, this appears to kill the model; after all, 
the extrapolation in the standard model leads to expect 
m T « 2777^. However, it is not possible to exclude that 



several effects modify this conclusion and avoid a flat 
contradiction with data (although we would in any case 
expect that mi, comes out as small as possible). In par- 
ticular, we note that with a suitable choice of phases the 
corrections order e can amount to a 10 % reduction, that 
the large Dirac Yukawa coupling can produce a 10 or 
20 % effect 44] , similar to what one can estimate for the 
change due to the full three flavor analysis. 

3) Last but not least there is an important relation be- 
tween the quark mixing V cb and the atmospheric mixing 
angle. Eq. (|18J1 shows that the main part of the (unphys- 
ical) up- and down-quark rotations are the same; thus, 
the quark mixing is found to be: 

\Vch\ = I Re£-icos20 A Im£| + 0(|e 2 |) (24) 

where £ = cos20a {^d — £[/)■ This equation demonstrates 
the successful coexistence of small and large mixing an- 
gles. In order for it to work quantitatively, 1008 20^1 
should be as large as possible, i.e. 9a should be as far 
as possible from the maximal value 45°. Strictly speak- 
ing, even this would not be sufficient if this prediction 
is taken at its face value. However the neglected effects 
from the first generation and the loops prevent us from 
sentencing this prediction and this model. 

The analysis we presented above could be in princi- 
ple changed by the two loop corrections in the Yukawa 
sector 41]. In the model I this consists only in the renor- 
malization of the original couplings, while in model II it 
could generate an effective coupling of the fermions to a 
one-index object (10fp) such as for example 



210 H 126 H 45 H 120 H , s 

16 F ^ — -16 F , 16 F — a r 16 F , ... (25) 



M 



M 



Such terms, even if present, are negligible in the present 
study of the two generation case, but they should be 
taken into account in the full analysis of the three gener- 
ations. 



III. SYMMETRY BREAKING PATTERNS AND 
NEUTRINO MASS 

As argued already in the introduction, SO(10) GUT 
works perfectly well without invoking supersymmetry. It 
is true that supersymmetry leads naturally to the uni- 
fication of gauge couplings, but the same effect can be 
equally achieved with left-right symmetry as an interme- 
diate scale. This is precisely what happens in SO (10). 
In the over-constrained models discussed in this paper, 
the Dirac neutrino Yukawa couplings are not arbitrary. 
Thus one must make sure that the pattern of intermedi- 
ate mass scale is consistent with a see-saw mechanism for 
neutrino masses. More precisely, the B — L-breaking scale 
responsible for right-handed neutrino masses cannot be 
too low. On the other hand, this scale, strictly speak- 
ing, cannot be predicted by the renormalization group 
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study of the unification constraints. The problem is that 
the right-handed neutrinos and the Higgs scalar responsi- 
ble for B — L breaking are Standard Model singlets, and 
thus have almost no impact (zero impact at one-loop) 
on the running. Fortunately, we know that the B — L 
breaking scale must be below SU(5) breaking, since the 
couplings do not unify in the Standard Model. Better 
to say, Mb-l < Mr, the scale of SU(2)# breaking, and 
hence one must make sure that Mr is large enough. This, 
together with proton decay constraints, will allow us to 
select between a large number of possible patterns of sym- 
metry breaking. 

Our task is simplified by the exhaustive study of sym- 
metry breaking in the literature, in particular the careful 
two-loop level calculations of Ref. [lflj |. Recall, though, 
that the (2,2, 15) field must lie below the GUT scale as 
discussed in sect. Ill ATI and although its impact on the 
running is very tiny, it must be included. 

The lower limit on Mr stems from the heaviest neu- 
trino mass 



M h 



(26) 



which gives Mr > 10 13 GeV or so. One can now turn 
to the useful table of Ref. [lfj, where the most general 
patterns of SO (10) symmetry breaking with two interme- 
diate scales consistent with proton decay limits are pre- 
sented. (Notice that the models with subscript 'b' in the 
table utilize 16h in place of 126h to break the SU(2)ij 
symmetry, and thus are not relevant for our discussion.) 

The above limit on Mr immediately rules out 
a number of the remaining possibilities; the most 
promising candidates are those with an intermediate 
SU(2) L xSU(2) i? x SU(4) c xP symmetry breaking scale 
(that is, PS group with unbroken parity). This is the 
case in which the breaking at the large scale is achieved 
by a Pati-Salam parity even singlet, for example con- 
tained in 54h- In the searching for a realistic symmetry 
breaking pattern one does not need to stick to the global 
minimum of the potential as in [Til ], but a local meta- 
stable minimum with a long enough lifetime will do the 
job as well. It has to be stressed however, that a big 
uncertainty is implicit in all models with complicated or 
unspecified Higgs sector, due to possibly large and un- 
controlled threshold corrections [45j . 

In any case, the nature of the GUT Higgs and the pat- 
tern of symmetry breaking will also enter into the fitting 
of fermion masses, since they determine the decomposi- 
tion of the light (fine-tuned) Higgs doublet (e.g., they 
provide relations among the parameters c\, C2, C3, C4 in 
eq. H12I) ). This point is often overlooked but it is essential 
in the final test of the theory. At this point, for us it is re- 
assuring that both the pattern of symmetry breaking and 
the nature of Yukawa interactions allow for a possibly re- 
alistic, predictive minimal model of non-supersymmetric 
SO(10). 



IV. SUMMARY AND OUTLOOK 

The recent years have witnessed an in-depth study of 
supersymmetric SO(10) grand unification based on the 
renormalizable see-saw mechanism. What has emerged 
is a possibly realistic picture for the unification of matter 
and forces with a predictive pattern for neutrino masses 
and mixings. The crucial point is that the SO (10) sym- 
metry may be sufficient by itself, without the need for any 
additional physics. While the theory has a number of ap- 
pealing features, it may be killed by its main ingredient: 
there may not be low-energy supersymmetry. It may be 
partially or completely broken. A nice example of partial 
breaking is the so-called split supersymmetry with light 
higgsinos and gauginos and heavy scalars. This picture 
allows for the interesting possibility of a radiative see- 
saw mechanism for neutrino masses, and another simple 
predictive version of the SO(10) theory. 

Since we know nothing about the existence of super- 
symmetry or the nature of its breaking, it is mandatory 
to study the non-supersymmetric version, as a part of the 
search for the SO(10) GUT. This was the scope of our pa- 
per. We have identified two potentially realistic, predic- 
tive Yukawa structures for the case of the renormalizable 
see-saw mechanism, based on a 126h- This choice is 
motivated by the fact that the alternative radiative see- 
saw seems to favor split supersymmetry [25]. We have 
focused on the renormalizable version simply in order to 
be predictive, without invoking unknown physics. 

The models require adding 10h or 120h fields. The 
latter is particularly interesting, due to the small num- 
ber of Yukawa couplings. Both models seem to require 
adding U(l )pq . While this may be appealing since it pro- 
vides the axion as a dark matter candidate, it is against 
the spirit of sticking to pure grand unification. 

A number of issues must be addressed in order to con- 
struct a fully realistic theory. The first task, as we re- 
peatedly argued, is a complete three generation numeri- 
cal study. This also includes the construction of the min- 
imal GUT Higgs sector and the study of its impact on the 
fermion masses and mixings. For a successful model, if 
any, one must study in turn the proton decay predictions, 
and in particular, the branching ratios that are calcula- 
ble in the over-constrained theories discussed here. This 
is a less urgent (but equally important) task simply due 
to the lack of experimental data. Beside proton decay, 
the other generic feature of grand unification is the exis- 
tence of magnetic monopoles which brings along the so- 
called monopole problem due to the over-production in 
the early universe. While there is always the possibility 
of the inflation solution, it is worth recalling that in non 
supersymmetric theories there are other interesting ways 
out of this impasse. These are for example the symmetry 
non-restoration at high temperature [46j and the possi- 
bility that unstable domain walls sweep the monopole 
away |47|. In principle, either of these alternative solu- 
tions can provide further constraints on the parameters 
of the theory. Yet another constraint comes from leptoge- 
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nesis, which finds its natural role in SO(10) with seesaw 
mechanism. 

Acknowledgments 

This work was completed during the 2005 Gran Sasso 
Summer Institute, and we wish to thank Z. Bcrezhiani 



for organizing an excellent meeting and for his warm hos- 
pitality. The work of G.S. was supported in part by Eu- 
ropean Commission under the RTN contract MRTN-CT- 
2004-503369; the work of B.B. by the Ministry of Educa- 
tion, Science and Sport of the Republic of Slovenia: the 
work of A.M. by CDCHT-ULA project No. C-1244-04- 
05-B. B.B. thanks Michal Malinsky and Miha Nemevsek 
for discussion. 



[1] H. Georgi, In Coral Gables 1979 Proceeding, Theory and 
experiments in high energy physics, New York 1975, 329 
and H. Fritzsch and P. Minkowski, Annals Phys. 93 
(1975) 193. 

[2] J. C. Pati and A. Salam, Phys. Rev. D 10 (1974) 275. 

[3] R. N. Mohapatra and J. C. Pati, Phys. Rev. D 11 (1975) 
2558; G. Senjanovic and R. N. Mohapatra, Phys. Rev. D 
12 (1975) 1502; G. Senjanovic, Nucl. Phys. B 153 (1979) 
334. 

[4] P. Minkowski, Phys. Lett. B 67 (1977) 421; T. Yanagida, 
proceedings of the Workshop on Unified Theories and 
Baryon Number in the Universe, Tsukuba, 1979, eds. 
A. Sawada, A. Sugamoto; S. Glashow, in Cargese 1979, 
Proceedings, Quarks and Leptons (1979); M. Gell-Mann, 
P. Ramond, R. Slansky, proceedings of the Supergrav- 
ity Stony Brook Workshop, New York, 1979, eds. P. 
Van Niewenhuizen, D. Freeman; R. Mohapatra, G. Sen- 
janovic, Phys.Rev.Lett. 44 (1980) 912. 

[5] G. Dvali and A. Vilenkin, Phys. Rev. D 70 

(2004) 063501 |arXiv:hep-th/0304043] ; G. Dvali, 
arXiv:hep-th/0410286 . 

[6] P. Langacker and M. x. Luo, Phys. Rev. D 44 (1991) 817. 
[7] S. Dimopoulos, S. Raby, F. Wilczek, Phys. Rev. D 24 
(1981) 1681. L.E. Ibanez, G.G. Ross, Phys. Lett. B 105 

(1981) 439. M.B. Einhorn, D.R. Jones, Nucl. Phys. B 196 

(1982) 475. W. Marciano, G. Senjanovic, Phys.Rev.D 25 
(1982) 3092. 

[8] C. S. Aulakh, B. Bajc, A. Melfo, G. Senjanovic 

and F. Vissani, Phy s. Lett. B 588, 196 (2004) 

|arXiv:hep-ph/0306242|. 

[9] G. Senjanovic, arXiv:hep-ph/0501244 
[10] N. G. Deshpand eTTT Keith an d P. B. P al, Phys. Rev. D 

47, 2892 (1993) [arXiv:hep-ph/9211232| . 
[11] F. Acampora, G. Amelino-Camelia, F. Buccella, O. Pi- 

sant i, L. Rosa and T. Tuz i, Nuovo Cim. A 108 (1995) 

375 |arXiv:hep -ph/9405332 . 
[12] R. D. Peccei and H. R. Quinn, Phys. Rev. Lett. 38 (1977) 

1440. 

[13] B. Bajc, A. Melfo, G. Senjanovic and F. Vissani, Phys. 

Rev. D 70 (2004) 035007 arXiv:hep-ph/0402122 . 
[14] C. S. Aula kh and A. Girdhar, Int. J. Mod. Phys. A 20 

(2005) 865 arXiv:hep-ph/0204097 . 

[15] T. Fukuyama, A. Ilakovac, T. Kikuchi, S. Meljanac 

and N. Okada, J. M ath. Phys. 46, 033505 (2005) 
arXiv:hep-ph/0405300 . 
[16] K. Matsuda, Y. Koide and T. Fukuyama, Phys. Rev. D 

64 (2001) 053015 |arXiv:hep- ph/0010026 . 
[17] K. Matsuda, Y. Koide, T. Fukuyama and H. Nishiura, 

Phys. Rev. D 65 (2002) 033008 [Erratum-ibid. 

D 65 (2002) 079904] [arXiv:hep-p h/0108202 ; 



T. Fukuy ama and N. Okada, JHEP 0211 (2002) 
011 [arXiv:hep -ph/0205066 . 
[18] H. S. Goh, R. N. Mohapatra and S. P. Ng , Phys. Lett. 
B 570 (2003) 215 arXiv:hep-ph/0303055 ; H. S. Goh, 
R. N. Mohapatra and S. P. Ng, Phys. Rev. D 68 (2003) 
115008 [arXiv:hep-ph /0308197 ; B. Dutta, Y. Mimura 
and R. N. M ohapatra, Phys. Rev. D 69 (2004) 115014 

arXiv:hep-ph/0402113 . 

[19] K~S~Babu and C. Macesanu, |arXiv:hep-ph /0505200 
[20] S. Bertolini and M. Malinsky, arXiv:hep-ph/0504241 
[21] B. Bajc, G. Senjanovic, F. Vissani, Phys. Rev. Lett. 90 

(2003) 051802 hep-ph/02 10207 , and hep-ph/0110310 
[22] B. Bajc, G. Senjanovic and F. Vissani, Phys. Rev. D 70 

(2004) 093002 arXiv:hep-ph/0402140 . 

[23] R. N. Mohapatra, Phys. Rev. D 34, 3457 (1986). 
A. Font, L. E. Ibanez and F. Quevedo, Phys. Lett. B228, 
79 (1989). S. P. Martin, Phys. Rev. D46, 2769 (1992). 

[24] C.S. Aulakh, K. Benakli, G. Senjanovic, Phys. Rev. 
Lett. 79 (1997) 2188. C. S. Aulakh, A. Melfo 
and G. Sen janovic, Phys. Rev. D 57, 4174 (1998); 
arXiv:hep-ph/9707256 . C. S. Aulakh, A. Melfo, 
A. Rasin and G. Senjanovic, Phys. Rev. D 58, 115007 
(1998) [arXiv:hep-ph/9712551 | . C. S. Aulakh, A. Melfo, 
A. Rasin and G. Senjanovic. Phys. Lett. B 459 (1999) 
557. |arXiv:hep-ph/9902409) . 

[25] B. Bajc and G. Senjanovic, arXiv:hep-ph/0507169 

[26] G. Lazarides, Q. Shan and C. Wetterich, Nucl. Phys. B 
181 (1981) 287. 

[27] R. N. Mohapatra and G. Senjanovic, Phys. Rev. D 23 
(1981) 165. 

[28] K. S. Babu and R. N. Mohapatra, P hys. Rev. Lett. 70 

(1993) 2845 |arXiv:hep-ph/9209 2T5l. L. Lavou ra, Phys. 

Rev. D48 (1993) 5440 arXiv:hep-ph/9306297 . 
[29] H. Arason, D. J. Castano, E. J. Piard and P. Ramond, 

Phys. Rev. D 47 (1993) 232 arXiv:hep-ph/9204225 . 
[30] S. Weinberg, Phys. Rev. Lett. 40 (1978) 223; F. Wilczek, 

Phys. Rev. Lett. 40 (1978) 279. 
[31] R. N. Mohapatra and G. Senjanovic, Z. Phys. C 17 

(1983) 53. 

[32] J. E. Kim, Phys. Rept. 150 (1987) 1. 
[33] R. Holman, G. Lazarides and Q. Shafi, Phys. Rev. D 27 
(1983) 995. 

[34] M. B. Wise, H. Georgi and S. L. Glashow, Phys. Rev. 

Lett. 47 (1981) 402. 
[35] J. R. Ellis and M. K. Gaillard, Nucl. Phys. B 150 (1979) 

141. 

[36] M. A. B. Beg and H. S. B. Tsao, Phys. Rev. Lett. 41 
(1978) 278. 

[37] R. N. Mohapatra and G. Senjanovic, Phys. Lett. B 79, 
283 (1978). 



8 



[38] H. Georgi, Hadronic J. 1 (1978) 155. 

[39] E. Witten, Phys. Lett. B 91 (1980) 81. 

[40] N. Arkani-Hamed and S. Dimopoulos, 

arXiv:hep-th/0405159 G. F. Giudice and A. Romanino, 

Nucl. Phys. B 699 (2004) 65 arXiv:hep-ph/0406088 . 

N. Arkani-Hamed, S. Dimopoulos, G. F. Giudice and 

A. Romanino, arXiv:hep-ph/0409232 
[41] B. Bajc and G. Senjano vic, Phys. Eett. B 610 (2005) 80 

arXiv:hep-ph/0411193 . 
[42] P. Nath and R. M. Syed, Phys. Lett. B 506 

(2001) 68 [Erratum -ibid. B 508 (2001) 216] 

arXiv:hep-ph/0103165 . 
[43] For a recent analysis and references, see A. Stru- 



mia and F. Vi ssani, Nucl. Phys. B 726 (2005) 294 
arXiv:hep-ph/0503246 . 
[44] F. Vissani and A. Y. Smirno v, Phys. Lett. B 341 (1994) 

173 arXiv:hep-ph/9405399 . 
[45] V. V. Dixit and M. Sher, Phys. Rev. D 40 (1989) 3765; 

C. S. Aulakh and A. Girdhar, Nucl. Phys. B 711 (2005) 

275 arXiv:hep-ph/0405074 . 
[46] G. R. Dvali, A. Melfo and G. Senjanovic, Phys. Rev. Lett. 

75, 4559 (1995) arXiv:hep-ph/9507230 . 
[47] G. R. Dvali, H. Liu and T. Vachaspati, Phys. Rev. Lett. 

80, 2281 (1998) arXiv:hep-ph/9710301 . 



